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Abstract 


Data  fusion  is  suitable  for  a  broad  range  of  decision  support  applications.  To  cope  with 
a  larger  class  of  problems  and  contexts,  data  fusion  gains  to  be  adaptive.  Adaptation  in 
data  fusion  corresponds  to  Level  4  of  the  JDL  model,  also  referred  to  as  process  refine¬ 
ment.  The  Decision  Support  Systems  Section  (DSS)  at  Defence  Research  &  Development 
Canada  (DRDC)  -  Valcartier  has  initiated  research  activities  aiming  at  developing  and 
demonstrating  advanced  concepts  of  adaptive  data  fusion  that  could  apply  to  the  current 
Halifax  and  Iroquois  Class  Command  &  Control  Systems  (CCS),  as  well  as  their  possi¬ 
ble  future  upgrades,  in  order  to  improve  their  performance  against  the  predicted  future 
threats.  This  document  gives  a  brief  description  of  the  adaptive  data  fusion  concepts.  It 
also  presents  a  new  Measure  Of  Effectiveness  (MOE)  that  serves  as  an  adaptation  trigger  in 
the  target-tracking  problem  in  maritime  Above  Water  Warfare  (AWW)  applications.  The 
proposed  MOE  uses  the  search  to  lock-on  time  of  the  Fire  Control  Radar  (FCR)  and  aims 
at  establishing  and  quantifying  the  effect  of  the  quality  of  the  Maritime  Tactical  Picture 
(MTP)  on  the  diminution  of  battle  space  size  and  reaction  time.  Besides  adaptation  of  the 
sensing  and  processing  operation,  this  MOE  allows  addressing  the  trade-off  between  the 
time  dedicated  to  the  tracking  with  surveillance  radars  versus  the  time  spent  in  search  and 
lock-on  with  FCR. 

Resume 


La  fusion  de  donnees  se  prete  a  une  vaste  gamme  duplications  pour  l’aide  a  la  decision. 
Afin  de  pouvoir  repondre  a  une  classe  etendue  de  problemes,  la  fusion  de  donnees  gagnerait  a 
etre  adaptative.  Ceci  definit  le  niveau  4  du  modele  JDL  de  fusion  de  donnees.  La  Section  des 
systemes  d’aide  a  la  decision  (SAD),  a  Recherche  et  developpement  pour  la  defense  Canada 
(RDDC)  -  Valcartier,  a  initie  des  activites  de  recherche  ayant  comme  principal  objectif  de 
demontrer  des  concepts  avances  de  fusion  adaptative  des  donnees.  Ces  concepts  pourraient 
s’appliquer  aux  systemes  actuels  de  commandement  et  controle  des  navires  canadiens  de 
classe  Halifax  et  Iroquois,  ainsi  qu’a  leurs  versions  futures.  Ce  document  fournit  une  breve 
description  des  concepts  de  fusion  adaptative  de  donnees.  Egalement,  une  nouvelle  mesure 
d’efficacite  est  definie  et  presentee.  Cette  mesure  servirait  de  detecteur  de  changement  dans 
la  situation  pouvant  activer  l’adaptation.  Cette  mesure  d’efficacite  emploie  le  temps  de 
recherche  du  radar  de  conduite  de  tir  et  vise  a  etablir  et  quantifier  l’effet  de  la  qualite  de 
l’image  tactique  maritime  sur  la  diminution  du  temps  de  reaction.  Cette  mesure  d’efficacite 
permet  egalement  d’aborder  le  probleme  de  la  recherche  de  compromis  entre  le  temps  alloue 
au  pistage  par  les  radars  de  surveillance  et  au  radar  de  conduite  de  tir  pour  la  recherche. 
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Executive  summary 


Impact  of  the  tactical  picture  quality  on  the  fire  control  radar 
search-lock-on  time 

F.  Rheaume,  A.  R.  Benaskeur ;  DRDC  Valcartier  TR  2003-284;  Defence  R&D 
Canada  -  Valcartier;  September  2006. 

Adaptive  data  fusion  is  aimed  at  developing  systems  that  can  cope  with  changes  in  their 
context  by  performing  structural  modifications.  This  constitutes  a  step  further  compared 
to  the  purely  passive  and  open  loop  data  fusion.  It  allows  obtaining  increased  quality  with 
equal  or  less  effort.  Adaptive  systems  can  be  seen  meeting  adaptation  requirements  that  are: 
goal  management,  performance  evaluation,  effectiveness  evaluation,  and  action  selection. 
These  requirements  should  enable  the  system  to  automatically  adapt  to  unforeseen  and 
changing  contexts,  thus  resulting  in  better  responses. 

This  report  describes  research  activities  conducted  at  Defence  Research  &  Development 
Canada  (DRDC)  -  Valcartier  that  aims  at  bringing  contributions  in  adaptive  data  fusion 
area  on  both  theoretical  and  experimental  sides.  The  objective  of  the  reported  work  is 
to  define  a  new  Measure  Of  Effectiveness  (MOE)  that  serves  for  the  adaptation  of  the 
target  tracking  operation  in  maritime  Above  Water  Warfare  (AWW)  context.  Measures  of 
Effectiveness  are  required  to  evaluate  the  system  behaviour  with  respect  to  high  level  and 
operational  objectives.  The  results  of  such  an  evaluation  set  the  low  level  objectives.  The 
proposed  MOE  uses  the  search  and  lock  on  time  of  the  Fire  Control  Radar  (FCR)  and  aims 
at  establishing  and  quantifying  the  effect  of  the  quality  of  the  Maritime  Tactical  Picture 
(MTP)  on  the  diminution  of  battle  space  size  and  reaction  time.  This  work  attempts  to 
demonstrate  the  benefits  of  adaptive  tracking  that  makes  use  of  reaction  time  as  a  context 
compared  to  the  classical  open  loop  tracking  mode. 

The  Fire  Control  Radar  is  involved  in  the  target  engagement  phase  of  the  maritime  AWW 
operations.  Once  a  decision  is  made  to  engage  a  given  target,  the  FCR  is  cued  with 
positional  information  from  the  surveillance-sensor-trackers  in  order  to  re-acquire,  track 
and  engage  it.  The  FCR  has  a  search  time  restricted  by  the  reaction  time  that  depends 
on  several  factors  such  as  the  ownship  weapons  properties,  Command  &  Control  System 
(CCS)  performance,  the  operator  skill/training,  and  by  the  attacking  target  characteristics. 
Besides  these  factors,  the  time  available  to  the  decision  is  affected  by  the  duration  of  search 
and  lock-on  operation.  This  time  the  FCR  will  take  to  re-acquire  the  target  is  directly 
dependent  on  the  size  of  uncertainty  volume  defined  in  the  tactical  picture.  A  poor  quality 
MTP  causes  the  FCR  to  search  in  a  large  size  volume.  Such  an  MTP  will  have  grave 
consequences  on  the  battle  space  size,  since  the  time  for  the  search  and  lock-on  depends,  in 
a  non-linear  manner,  on  the  volume  to  be  scanned  and  the  probability  of  detection  of  the 
FCR.  This  time  is  subtracted  from  the  total  reaction  time  necessary  to  the  ship’s  survival. 
This  search  to  lock-on  time  constitutes  the  new  Measure  Of  Effectiveness  proposed  in  this 
document.  A  relationship  between  the  track  quality  and  the  diminution  of  battle  space  size 
and  reaction  time  is  defined. 
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La  fusion  adaptative  de  donnees  vise  a  developper  des  systemes  capables  de  repondre 
adequatement  aux  changements  de  leur  contexte.  Ce  concept  constitue  une  evaluation  par 
rapport  a  la  fusion  passive  et  en  boucle  ouverte.  Les  systemes  adaptatifs  peuvent  etre  percus 
comme  repondant  aux  exigences  d’adaptation  que  sont  la  gestion  des  buts,  1’evaluation  de 
la  performance,  revaluation  de  l’efficacite  et  la  selection  des  suites  d’actions. 

Ce  rapport  decrit  les  activites  de  recherche,  menees  a  Recherche  et  developpement  pour  la 
defense  Canada  (RDDC)  -  Valcartier,  qui  visent  a  apporter  une  contribution  au  domaine 
de  la  fusion  adaptative,  tant  sur  le  plan  theorique  qu’experimental.  L’objectif  du  travail 
presente  dans  ce  document  est  de  proposer  une  nouvelle  mesure  d’efficacite  qui  peut  agir 
comme  declencheur  d’adaptation  dans  les  applications  de  pistage  de  cible.  Pour  les  systemes 
adaptatifs,  les  mesures  d’efficacite  sont  souvent  requises  pour  evaluer  les  performances  du 
systeme  par  rapport  aux  objectifs  de  haut  niveau.  Le  resultat  de  cette  evaluation  etablit  les 
objectifs  de  niveau  inferieur.  La  mesure  d’efficacite  proposee  utilise  le  temps  de  recherche 
du  radar  de  conduite  de  tir  et  vise  a  etablir  et  a  quantifier  la  relation  de  dependance  entre 
la  qualite  de  l’image  tactique  et  la  diminution  de  temps  de  reaction.  Le  travail  presente 
tente  de  demontrer  l’avantage  du  pistage  adaptatif  qui  utilise  le  temps  de  reaction  comme 
contexte,  compare  au  pistage  classique  en  boucle  ouverte  sans  retroaction. 

Dans  le  cadre  de  guerre  anti-aerienne  et  de  surface,  le  radar  de  conduite  de  tir  fait  partie 
de  la  phase  d’engagement  des  cibles.  Une  fois  la  decision  prise  pour  engager  une  cible,  ce 
radar  recoit  du  systeme  de  pistage  les  informations  cinematiques  de  la  cible.  II  entame  alors, 
sur  la  base  de  cette  information,  une  operation  de  recherche  afin  de  reacquerir  et  pister  la 
cible  a  son  tour  pour  fins  d’engagement.  Le  temps  de  recherche  du  radar  de  tir  est  limite 
par  le  temps  de  reaction  qui,  a  son  tour,  depend  de  plusieurs  facteurs.  II  faut  noter  que 
cette  dependance  est  reciproque.  Le  temps  de  reaction  disponible  au  decideur  est  egalement 
affecte  par  la  duree  de  la  recherche  par  le  radar  de  conduite  de  tir.  Le  temps  que  met 
le  radar  de  conduite  de  tir  pour  reacquerir  la  cible  depend  egalement  directement  de  la 
taille  du  volume  d’incertitude  que  definit  l’image  tactique.  Une  qualite  inferieure  de  l’image 
tactique  oblige  le  radar  de  conduite  de  tir  a  chercher  dans  un  large  volume.  Une  qualite 
de  l’image  tactique  peut  reduire  considerablement  le  temps  de  reaction,  puisque  l’operation 
de  recherche  depend  d’une  facon  non  lineaire  du  volume  a  balayer  et  de  la  probability  de 
detection  du  radar  de  conduite  de  tir.  Cette  duree  est  soustraite  du  temps  de  reaction 
total  disponible  au  decideur  et  necessaire  a  la  survie  du  bateau.  La  duree  de  la  phase  de 
recherche  du  radar  de  conduite  de  tir  constitue  la  nouvelle  mesure  d’efficacite  proposee  dans 
ce  document.  La  relation  entre  la  qualite  des  pistes  et  la  diminution  du  temps  de  reaction 
est  definie. 
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1  Introduction 


Modern  ship-borne  Command  and  Control  (C2)  Systems  should  make  increasing  use  of  data 
fusion  technology  and  tools.  By  reducing  uncertainty  in  the  existing  pieces  of  information 
and  providing  means  to  infer  about  the  missing  pieces,  data  fusion  supports  the  operators 
in  compiling  and  analyzing  the  tactical  picture  and  ultimately  improving  the  situation 
awareness  of  the  decision  makers. 

Data  fusion  has  often  been  portrayed  as  a  set  of  sequential  processes  that  combine  pieces  of 
information  in  a  purely  passive  and  open  loop  fashion.  Nevertheless,  to  cope  with  changing 
objectives,  environments  and  constraints,  and/or  to  achieve  a  higher  performance,  modern 
C2  systems  need  many  additional  functions.  The  most  essential  of  them  is  an  active  feedback 
or  adaptation  [1].  The  introduction  of  such  a  functionality  into  the  situation  analysis 
process  and  the  underlying  data  fusion  processing  defines  a  part  of  the  emerging  discipline 
of  process  refinement  (or  adaptive  data  fusion).  However,  even  though  it  has  witnessed  a 
growing  interest  during  the  last  few  years,  the  process  refinement,  that  represents  a  logical 
extension  of  the  primary  data  fusion  concepts,  still  lacks  mature  results  and  a  shared  precise 
definition. 

The  Decision  Support  Systems  (DSS)  Section,  at  Defence  Research  &  Development  Canada 
(DRDC)  -  Valcartier,  is  conducting  research  activities  aiming  at  defining,  developing  and 
demonstrating  advanced  concepts  of  data  fusion  adaptation  and  sensor  management  that 
could  apply  to  the  current  Command  &  Control  Systems  (CCS)  of  the  Halifax  and  Iroquois 
Class  ships  as  well  as  their  possible  future  upgrades. 

This  report  describes  research  activities  that  aims  at  bringing  contributions  on  both  the¬ 
oretical  and  experimental  sides.  The  objective  here  is  to  define  a  new  Measure  Of  Effec¬ 
tiveness  (MOE)  that  serves  for  the  adaptation  of  the  target  tracking  operation  in  maritime 
Above  Water  Warfare  (AWW)  context.  Note  that,  in  adaptive  systems,  Measures  of  Effec¬ 
tiveness  are  often  required  to  evaluate  the  system  behaviour  with  respect  to  high  level  and 
operational  objectives.  The  results  of  such  an  evaluation  set  the  low  level  objectives.  The 
proposed  MOE  uses  the  search  to  lock-on  time  of  the  Fire  Control  Radar  (FCR)  and  aims 
at  establishing  and  quantifying  the  effect  of  the  quality  of  the  Maritime  Tactical  Picture 
(MTP)  on  the  diminution  of  battle  space  size  and  reaction  time.  This  work  attempts  to 
demonstrate  the  benefits  of  adaptive  tracking  that  makes  use  of  reaction  time  as  a  context 
compared  to  the  classical  open-loop  tracking  mode. 

The  Fire  Control  Radar  is  involved  in  the  target  engagement  phase  of  the  maritime  AWW 
operations.  Once  a  decision  is  made  to  engage  a  given  target,  the  FCR  is  cued  with 
positional  information  from  the  surveillance-sensor-trackers  in  order  to  re-acquire,  track, 
and  engage  it.  The  FCR  has  a  search  time  restricted  by  the  reaction-time  that  depends 
on  several  factors  such  as  the  ownship  weapons  properties,  Command  &  Control  System 
(CCS)  performance,  the  operator  skill/training,  and  by  the  attacking  target  characteristics. 
Besides  these  factors,  the  time  available  to  the  decision  is  affected  by  the  duration  of  search 
and  lock-on  operation.  This  time  the  FCR  will  take  to  re-acquire  the  target  is  directly 
dependent  on  the  size  of  uncertainty  volume  defined  in  the  tactical  picture.  A  poor  quality 
MTP  causes  the  FCR  to  search  in  a  large  size  volume.  Such  an  MTP  will  have  grave 
consequences  on  the  battle  space  size,  since  the  time  for  the  search  and  lock-on  depends,  in 
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a  non-linear  manner,  on  the  volume  to  be  scanned  and  the  probability  of  detection  of  the 
FCR.  This  time  is  subtracted  from  the  total  reaction  time  necessary  to  the  ship  survival. 

This  search  to  lock-on  time  constitutes  the  new  Measure  Of  Effectiveness  proposed  in  this 
document.  A  relationship  between  the  track  quality  and  the  diminution  of  battle  space  size 
and  reaction  time  is  defined. 

The  proposed  MOE  can  be  used  as  an  adaptation  trigger  in  many  problems  in  which  the 
fire  control  radar  is  involved.  A  simple  adaptation  problem  consists  in  turning  the  tracking 
operation  on/off  depending  on  the  time  available  and  the  required  quality  of  the  tracks.  A 
more  general  problem  is  to  adapt  the  tracking  by 

1.  determination  of  the  tracking  duration  before  cueing  the  FCR 

2.  determination  of  the  tracking  duration  and  nature  (adapt  processing  and  resources) 
to  meet  the  requirements. 

Finally,  the  problem  of  establishing  relationship  between  the  track  quality  and  the  diminu¬ 
tion  of  the  battle  space  size  and  the  reaction  time  is  considered.  This  allows  defining  time 
constraints  for  the  search  and  lock-on  operation,  given  the  tactical  situation.  A  maximum 
bound  on  the  search  to  lock-on  time  is  determined  knowing  the  reaction  time  available  to 
the  decision  maker. 

This  report  is  organized  as  follows.  Section  2  introduces  concepts  and  definitions  of  adaptive 
systems.  This  section  is  intended  to  situate  the  MOE  computation  task  within  the  adap¬ 
tation  loop.  Section  3  presents  the  fire  control  problem,  within  the  whole  detect-to-engage 
sequence.  The  search  and  lock-on  problem  is  also  described.  In  Section  4,  the  time  estima¬ 
tion  process  is  presented.  The  memorandum  is  concluded  in  Section  5.  Appendix  A  gives 
the  state  prediction  equations,  while  in  Appendix  B  the  coordinate  conversion  equations 
are  given. 
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2  Adaptive  systems 


A  system  is  called  adaptive  when  it  can  analyze  its  own  performance  and  dynamically 
reconfigures  itself  to  compensate  for  changes  in  its  context.  The  changes  might  be  from 
the  environment  it  evolves  in,  its  objectives  and/or  its  needs.  Instead  of  being  developed 
for  a  specific  situation,  adaptive  systems  are  therefore  able  to  handle  a  family  of  situations 
defined  by  a  set  of  structural  constraints  on  their  context.  Designers  cannot  foresee  all 
the  circumstances  the  systems  will  be  used  in.  Therefore,  adaptation  mechanisms  aim  at 
making  those  systems  automatically  adapt  to  unforeseen  and  changing  contexts,  by  allowing 
for 

1.  changing  the  system  behaviour  and/or  performance  without  the  need  to  entirely  re¬ 
design  it, 

2.  doing  more  on  behalf  of  the  user  without  a  constant  interaction,  and 

3.  yielding  higher  performance  or  lower  costs,  or  both  whenever  possible. 

An  adaptive  system  perceives  what  is  going  on  around  itself1  and  then  uses  this  knowledge 
to  produce  appropriate  actions  to  achieve  its  goals.  The  system  must  take  into  account 
the  dynamic  nature  of  the  environment  it  interacts  with  to  reach  its  goals  throughout 
a  wide  range  of  situations  with  the  desired  level  of  performance2.  Based  on  the  newly 
available/inferred  information,  the  adaptive  system  needs  to  modify  its  operation  mode  to 
cope  with  different  performance  criteria,  variable  contextual  conditions  and  changing  status 
of  resources. 

2.1  Adaptation  requirements 

Adaptive  systems  are  aimed  at  producing  applications  that  can  readily  adapt  to  changing 
operating  environment  and  user’s  needs  and  desires.  Adaptive  systems  need  to  explicitly 
represent  the  actions  that  can  be  taken  and  the  goals  that  the  user  is  trying  to  achieve. 
This  makes  possible  for  the  user  to  change  goals  without  the  need  to  redesign  the  system. 
Also,  the  system  does  more  on  behalf  of  the  user  without  constant  interaction. 

An  adaptive  system  must  therefore  be  able  to  measure  and  detect  changes  in  its  perfor¬ 
mance  and  respond  to  these  changes  by  performing  structural  changes  either  parametric 
(configuration),  function/reasoning-path  level  (processing)  or  at  the  process  level  (resource- 
bounded  reasoning).  See  [2]  for  a  more  detailed  discussion  of  adaptation  requirements  and 
mechanisms. 

Therefore,  policies  that  can  dynamically  perform  these  structural  changes  in  response  to 
a  changing  context  are  required.  There  are  major  requirements  an  adaptive  system  must 
meet  to  adequately  operate  in  its  environment.  These  requirements  are  given  in  Figure  1 
and  discussed  in  the  followings. 

To  obtain  some  idea  of  the  state  of  its  environment,  namely,  its  context. 

2which  may  be  constant  or  variable. 
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Constraints 


Figure  1:  Adaptation  Requirements 


2.1.1  Goal  management 

Handling  the  adaptation  problem  presumes  that  the  required  performance  of  the  system 
can  be  specified  quantitatively  to  define  the  adaptation  goals/objectives.  Metrics  can  then 
be  calculated  and  used  to  evaluate  the  system’s  performance  (i.e.,  the  deviation  from  the 
desired  behaviour).  Another  important  requirement  is  to  keep  the  adaptation  effort  (time 
and  resources  consumption)  minimal. 

Generally,  specifications  can  be  divided  into  two  categories:  performance  specifications 
and  robustness  specifications.  Both  need  to  be  quantified  to  specify  the  adaptation  goals. 
Although  the  boundaries  between  the  two  can  be  fuzzy,  the  performance  specifications 
describe  the  desired  response  of  the  nominal  system,  i.e.,  in  the  absence  of  uncertainty. 
Nevertheless,  for  a  given  system,  the  nominal  performance  alone  may  not  suffice.  The 
system  must  keep  providing  acceptable  performance  even  if  the  context  is  different  (in 
some  measured  way)  from  the  intended  one.  Therefore,  robustness  specifications  are 
required  to  limit  the  degradation  of  the  performance  in  the  presence  of  uncertainty  that 
may  come  in  three  basic  versions: 

1.  uncertainties  in  the  data  (or  non-parametric  uncertainties), 

2.  uncertainties  in  the  system  models  (also  known  as  parametric  uncertainties),  and 

3.  constraints  on  the  admissible  actions  (the  span  of  the  adaptation  space). 

It  is  worth  mentioning  that  goal  management  is  not  required  in  all  adaptive  systems  since 
adaptation  objectives  are  often  static.  The  desired  level  of  performance,  or  reference  tra¬ 
jectory  that  the  adaptive  system  tries  to  achieve,  is  specified  beforehand  and  remains  un¬ 
changed.  The  difference  between  this  reference  trajectory  and  the  measured  level  of  per¬ 
formance  provides  a  good  measure  of  the  system  actual  behaviour  with  respect  the  desired 
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one.  These  metrics  are  used  as  an  action  selection  (see  Subsection  2.1.4)  basis  to  reduce 
any  observable  discrepancy. 

2.1 .2  Performance  evaluation 

Given  the  knowledge  of  the  current  state  of  the  environment  and  the  adaptation  objective 
(i.e.,  the  context),  a  very  important  aspect  of  adaptive  systems  concerns  the  performance 
evaluation  issue.  This  consists  in  finding  a  Measure  of  Performance  (MOP)  that,  once 
optimized,  leads  to  the  most  desirable  outcome.  The  selected  MOP  serves  as  an  action 
selection  basis  (or  utility  function  in  decision  theory  terminology).  It  is  required  to  grade 
the  benefits  from  the  different  possible  actions  so  that  an  optimal  solution  can  be  chosen. 
For  instance,  in  the  target-tracking  context,  the  information  update  paradigm  may  lead 
to  an  intuitive  method  of  addressing  the  MOP  selection  problem.  The  adaptive  system 
should  therefore  be  able  to  monitor  its  behaviour  and  detect  deviations  from  commitments 
or  new  opportunities.  It  should  also  be  able  to  accept  new  needs  from  external  sources  and 
evaluate  deviations  with  respect  to  current  commitments. 

2.1 .3  Effectiveness  evaluation 

While  the  MOPs  measure  deviations  from  commitments  on  low  (system)  level  objectives, 
a  set  of  Measure  of  Effectiveness  (MOE)  is  required  to  evaluate  the  system  behaviour  with 
respect  to  high  (operational)  level  objectives.  The  results  of  such  evaluation  define  the  low 
level  objectives  (for  the  inner  loop).  The  MOEs  could  depend  on  one  or  several  MOPs. 

The  search  to  lock-on  time  we  are  interested  in  this  document,  which  affects  the  battle  space 
size  and  reaction  time,  gives  a  good  example  of  a  possible  MOE  in  naval  AWW  context. 

2.1.4  Action  selection 

Given  the  goal  specifications,  the  environment  changes  and  the  performance  and/or  effec¬ 
tiveness  measures,  the  adaptation  problem  amounts  to  selecting  the  appropriate  course  of 
actions.  The  objective  of  these  actions  is  to  maintain  the  discrepancy  between  the  measured, 
or  estimated,  current  performance  and  the  desired  one  within  some  fixed  bounds.  This  dis¬ 
crepancy  (or  error)  measures  how  the  system  being  considered  behaves  with  respect  to  the 
specified  level  of  performance.  Therefore,  adaptation  mechanisms  should  be  able  to  reason 
and  make  commitments  on  the  goal/environment  changes.  Different  techniques  exist  for 
action  selection.  The  action  selection  defined  above  is  the  essence  of  the  control  design  that 
consists  in  selecting  the  input  (actions)  to  a  process  so  that  its  output  follows  the  desired 
specifications,  while  not  requiring  too  much  control  effort  (time/resource  consumption). 

When  no  control  is  used,  the  process  input  is  selected  at  the  time  of  design  and  remains  the 
same  throughout  the  lifetime  of  the  process.  When  control  is  used,  there  are  two  sources  for 
its  activation.  These  are  the  goal  modifications  and  relevant  environmental  changes.  Each 
one  corresponds  to  a  different  control  strategy.  Feedforward  control  is  concerned  with  the 
goal  changes,  while  feedback  control  handles  changes  in  the  environment.  The  interactions 
between  these  two  components  form  a  closely  integrated  loop  for  systems  to  adapt  to  its 
challenges.  In  the  Artificial  Intelligence  (AI)  planning  context,  these  respectively  correspond 
to  the  purely  deliberative  planning  model  and  the  purely  reactive  planning  model. 
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Once  a  course  of  actions  is  selected,  configuration  modification  functionality  must  offer 
to  the  adaptive  system  the  means  for  its  execution.  The  system  should  be  able  to  insert 
the  change  with  minimal  disruption  to  existing  behaviours.  Two  properties  are  of  prime 
importance  when  executing  actions.  These  are: 

1.  Stability  —  It  is  the  most  important  parameter  when  dealing  with  closed-loop  sys¬ 
tems.  A  system  is  said  stable  if  its  response  to  a  bounded  input  is  itself  bounded  by 
a  desirable  range. 

2.  Response  time  —  The  response  (transient)  time  measures  the  transient  response  of 
the  system.  That  is  the  time  the  system  takes  to  transit  from  one  stable  state  to 
another. 
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3  Problem  statement 


In  order  to  engage  a  given  target,  the  trackers  associated  with  the  search  and  surveillance 
sensors  (or  tracks  that  result  from  their  fusion)  cue  the  FCR.  The  accuracy  of  the  infor¬ 
mation  provided  to  the  FCR  determines  the  volume  it  must  scan  in  order  to  re-acquire  the 
target  and  lock  on  it  (see  Figure  2).  The  time  the  FCR  will  take  to  re-acquire  the  target  is 
directly  dependent  on  the  size  of  that  volume.  It  may  also  happen  that  the  FCR  will  never 
re-acquire  the  target.  A  poor  quality  MTP,  i.e.,  with  large  uncertainty  regions,  causes  the 
FCR  to  search  in  a  large  size  volume.  Such  an  MTP  will  have  grave  consequences  on  the 
battle  space  size.  The  time  for  the  search  and  lock-on  depends,  in  a  non-linear  manner, 
on  the  volume  to  be  scanned  and  the  probability  of  detection  of  the  FCR.  This  time  is 
subtracted  from  the  total  reaction  time  available  to  the  decision  maker  and  necessary  to 
the  ship  survival. 


Figure  2:  FCR  Search  Volume 


3.1  Detect-to-engage  sequence 

The  Fire  Control  Radar  is  involved  in  the  target  engagement3  phase  of  the  maritime  Above 
Water  Warfare  (AWW)  operations.  Prior  to  this  engagement  phase  there  is  a  threat  eval¬ 
uation  process  that  starts  with  detection/tracking  and  ends  with  the  engagement  deci- 
sion/FCR  tasking. 

Upon  detection  by  the  2D  search  and  surveillance  radars,  the  contact  information  is  provided 
to  the  tracking  system  that  maintains  a  more  accurate  estimation  of  the  target  position. 
Target  identity  inference  process  is  executed  concurrently. 


In  the  sequel,  only  the  target  position  will  be  considered.  It  is  given  by  state  estimate 

-I  T 


X  = 


x,y,x,y 


(i) 


3It  may  happen  that  the  FCR  be  used  in  the  Situation  Analysis  Process  to  obtain  3D  information  of  a 
given  target  or  to  make  the  target  uncover  its  actual  indent. 
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and  the  related  covariance  matrix  P.  This  information  is  used  by  subsequent  processes  to 
evaluate  the  threat  level  of  the  different  targets  present  in  the  volume  of  interest. 

Once  an  engagement  decision  is  made,  the  kinematical  information  of  the  target  to  proceed 
is  used  to  cue  the  FCR,  i.e.,  to  delimit  its  search  region.  The  FCR  initiates  its  search 
within  an  area  defined  mostly  by  the  2D  bearing- elevation  information  obtained  from  the 
Cartesian  position  (x,P)  provided  by  the  tracking  system.  Following  a  specific  pattern, 
the  FCR  will  scan  this  region  of  the  volume  of  interest  until  it  detects  and  locks  on  the 
target  for  which  a  track  is  then  maintained.  The  target  course  and  speed  contained  in  this 
FCR  track  is  then  used  to  compute  a  Predicted  Intercept  Point  (PIP)  inside  the  weapon 
engagement  envelope. 

The  problem  investigated  here  is  to  establish  quantitatively  the  impact  of  the  track  quality 
(defined  by  P )  on  the  search  to  lock-on  time  of  the  FCR. 

3.2  FCR  characteristics 

During  the  target  detection  phase,  the  search  and  surveillance  radars  scan  the  volume  of 
interest  for  contacts.  For  Canadian  Patrol  Frigates  (CPFs),  search  and  surveillance  radars 
(i.e.,  SPS-49  and  SG-150)  work  in  a  2D  space  and  therefore  provide  bearing  (/?)  and  range 
(R)  information,  as  shown  in  Figure  3.  No  elevation  is  provided  by  2D  radar. 

Note  that  the  bearing  illustrated  in  Figure  3  is  a  relative  bearing  rather  than  the  true 
bearing.  The  relative  bearing  is  the  angle  between  the  centre-line  of  the  ship  and  a  line 
pointed  at  the  target  measured  clockwise  from  the  bow.  The  true  bearing  is  the  angle 
between  the  true  north  and  a  line  pointed  at  the  target  measured  clockwise  from  the  true 
north.  Figure  3  also  shows  the  Cartesian  coordinate  system  with  the  X  and  Y  axes  positions. 


X 

Figure  3:  Search  and  surveillance  radar  2D  scanning 


Besides  the  fact  that  it  can  provide  the  elevation  (and  therefore  operates  in  3D),  the  FCR 
(such  as  the  CPF’s  STIR)  differs  from  the  search  and  surveillance  radars  by  its  circular  and 
much  narrower  beam.  This  results  in  increased  bearing  and  elevation  precision.  It  also  has 
higher  transmitting  frequency,  so  that  more  accurate  range  information  is  obtained. 
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Figure  4  shows  the  coordinate  system  for  the  FCR  along  with  the  range  (R),  the  relative 
bearing  (/?)  and  the  elevation  (e)  measurements.  The  X  and  Y  axes  positioning  is  identical 
to  2D  radars  where  the  Y  axis  points  at  the  bow  of  the  ship.  Because  most  of  the  tracking 
systems,  either  they  track  in  2  or  3  dimensions,  work  in  Cartesian  coordinates,  the  spherical 
measures  are  often  converted  into  Cartesian  ones  and  vice  versa.  This  allows  data  to  comply 
with  the  FCR  and  the  tracking  systems  needs  (see  Appendix  B  for  detailed  conversion 
equations). 


Figure  4:  Fire-control  radar  (FCR)  scanning  and  coordinate  system 

3.2.1  FCR  parameters 

The  FCR  is  characterized  by  its  beam  dimensions.  This  is  defined  by  its  width  bw  (or 
diameter  bw)  and  by  its  displacement  speed  b£  in  elevation  and  bp  in  bearing.  The  FCR 
has  a  search  path  that  can  follow  different  patterns.  This  search  path  must  be  taken  into 
account  in  the  estimation  of  the  search  to  lock-on  time. 


Characteristics 

Symbol 

Beam  diameter 

bw 

Beam  speed 

be,  b0 

Conditional  detection  probability 

Pd\l(x,v) 

Search  path 

~ 

Table  1:  FCR  parameters 


3.3  Search  to  lock-on  time 

Once  a  decision  is  made  to  engage  the  target,  the  FCR  is  cued  with  positional  information 
from  the  surveillance-sensor-trackers  in  order  to  re- acquire,  track,  and  engage  it.  The  FCR 
has  a  search  time  restricted  by  the  reaction  time  that  depends  on  several  factors.  This  might 
be  the  ownship  weapons  properties,  Command  &  Control  System  (CCS)  performance,  the 
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operator  skill/training,  and/or  the  attacking  target  characteristics.  Given  its  high-risk 
consequences,  the  engagement  phase  requires  more  precise  information  than  the  threat 
evaluation  phase.  This  is  why  the  3D  Fire  Control  Radar  takes  over  the  less  accurate  2D 
search  and  surveillance  radars.  To  provide  such  accurate  information,  the  FCR  need  to 
track  the  target  by  itself.  This  operation  cannot  be  performed  without  a  re-acquisition  that 
involves  a  search  and  lock-on  phase.  This  phase  starts  once  the  FCR  is  cued  by  the  tracking 
system  and  begins  its  scan.  It  ends  when  the  FCR  locks  on  the  target.  The  estimation  of 
the  time  (dsi)  the  FCR  will  take  to  perform  this  re- acquisition  task  defines  the  core  of  the 
addressed  problem. 

3.3.1  Applications 

The  new  investigated  MOE  aims  at  establishing  a  quantitative  relationship  between  the 
track  quality  and  the  search-lock-on  time  for  the  FCR.  Once  this  MOE  is  fed  back  (see 
Figure  5)  to  the  goal  manager,  different  adaptation  problems  that  exploit  this  contextual 
information  can  be  defined.  The  output  of  the  goal  manager  will  provide  the  adaptation 
triggers  for  the  inner-loop  that  acts  on  the  tracking  algorithms  and  the  resource  allocation  in 
order  to  achieve  the  specified  objectives.  The  inner-loop  measures  the  current  performance 
of  the  tracking  systems,  compares  this  performance  with  the  specified  goal,  and  selects 
appropriate  actions  to  reduce/eliminate  any  observed  discrepancy. 


Figure  5:  FCR  tpcR  based  tracking  adaptation  loop 


Examples  of  adaptation  problems  currently  investigated  and  that  make  use  of  search  to 
lock-on  time  as  a  trigger4,  are  defined  below. 

Problem  1:  Given  the  available  reaction  time,  the  quality  of  the  sensor  reports, 
the  delay  created  by  the  tracking  processing,  and  the  above  described  MOE,  one  can 
determine  whether  it  is  more  interesting  to  cue  the  FCR  directly  with  the  sensor 
reports  or  to  spend  more  time  to  improve  the  quality  of  the  information  (through  the 

4See  [2],  for  definitions. 
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tracking  operation)  before  cueing  the  FCR.  This  adaptation  problem  boils  down  to 
put  the  tracking  operation  on/off  depending  on  the  context. 

Problem  2:  Another  more  general  and  more  interesting  problem  consists  in  deter¬ 
mining,  depending  on  the  same  type  information  as  for  the  first  problem,  the  amount 
of  time  the  tracking  module  should  spend  in  improving  the  quality  of  the  information 
before  cueing  the  FCR.  This  problem  goes  beyond  just  switching  on/off  the  tracking 
function. 


Problem  3:  The  two  problems  described  above  assume  fixed  configurations.  The 
adaptation  in  the  first  one  boils  down  to  decide  whether  it  is  more  interesting  to  track 
or  not.  The  second  consists  in  calculating  the  best  duration  for  the  tracking  operation. 
The  problem  described  here  goes  further  and  assumes  redundancy  in  both  sensing  and 
processing  capabilities.  Hence,  besides  the  FCR,  the  ship  has  two  different  search  and 
surveillance  sensors.  One  sensor  is  assumed  to  provide  more  accurate  contacts  than 
the  other.  There  is  also  a  local  tracker  associated  with  each  sensor.  These  trackers  use 
a  set  of  different  tracking  algorithms,  each  providing  a  different  tracking  performance. 
Now  given  the  available  reaction  time,  the  quality  of  the  sensor  reports,  the  delay 
created  by  the  tracking  operation,  and  the  search  to  lock-on  time  MOE,  one  can 
dynamically  decide  to 

(a)  Switch  on/off  the  tracking  (problem  1) 

(b)  Use  one  sensor  for  tracking  with  variable  duration  (problem  2). 

(c)  Use  the  best  of  the  two  trackers  or  fuse  them 

(d)  Select  among  the  set  of  available  algorithms  for  tracking/ fusion  the  most  appro¬ 
priate  given  the  context. 

(e)  Combine  the  above  listed  problems  in  different  ways. 

Problem  4:  This  problem  aims  at  establishing  time  constraints  for  the  search  and 
lock-on  operation.  This  allows,  given  a  tactical  situation,  calculating  the  maximum 
affordable  time  for  search  and  lock-on  phase.  Given  the  reaction  time  available  to 
the  decision  maker,  the  latter  can  determine  the  maximum  search  to  lock-on  time  he 
can  afford.  Since  this  time  depends  on  the  quality  of  the  tactical  picture,  the  time 
requirements  can  be  re-expressed  as  quality  requirements. 

It  is  assumed  that  the  ownship  defensive  weapons  W  (such  as  SAM  missiles)  have 
an  effective  zone  between  a  minimum  range  Rw_  and  a  maximum  range  Rw+.  The 
coverage  zone  of  the  FCR  that  guides  those  weapons  is  defined  by  a  minimum  range 
Rf_  and  a  maximum  range  Rf+  (see  Figure  6). 

Let  Rw  be  the  defensive  weapon  speed  and  tw_  the  time  this  weapon  takes  before 
entering  its  effective  zone  (at  range  Rw_).  Then 

.  Rw 


Given  tw_  and  Rw_ ,  there  is  a  minimum  distance  Rt_  from  the  ship  (see  Equation  3) 
where  the  threat  R  must  be  detected  by  the  FCR  and  engaged  by  W . 


(2) 


Rf, _  —  Rw_  Rwtn 


(3) 
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Ownship 


Rf 


Defensive  weapon  effective  zone 

r  ^ 


R 
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Ru 


R, 


Fire  control  radar  coverage  zone 


T 


Figure  6:  Scenario  characteristics 


where  Rw  <  0.  If  the  defensive  weapon  W  is  launched  while  the  threat  T  has  already 
passed  Rt_,  it  will  be  too  late.  The  defensive  weapon  will  not  reach  its  target,  since 
the  latter  will  be  outside  its  effective  zone.  Therefore,  the  weapon  W  must  be  fired 
while  the  threat  R  is  beyond  Rt_.  The  time  available  to  the  FCR  (tfa)  is  given  by 


tfa  ~  tw- 


Rtc 

Rt 


(4) 


where  Rt  <  0  and  Rtc  is  the  range  of  the  threat  when  the  FCR  is  cued.  The  validity 
interval  of  this  time  is  given  by 


Rf+  >  Rtc  >  max{i?j_ ,  Rw_  +  dsiRt} 


(5) 


where  dsi  is  the  duration  of  the  search  and  lock-on  operation  of  the  FCR.  The  duration 
is  limited  by  the  available  time  tfa  given  by  Equation  4.  Obviously,  the  available  time 
is  maximum  when  the  FCR  is  cued  at  Rf+ 


tfa  ~  tw- 


Rf+  Rw— 
Rt 


(6) 


The  search  and  lock-on  duration  dsi  represents  the  core  problem  discussed  in  this 
document. 


One  can  easily  see  that  several  adaptation  problems,  that  exploit  the  high  level  contextual 
information  provided  by  the  search  and  lock  on  time,  can  be  defined  and  used  to  improve 
the  quality  of  the  response  of  the  ownship  to  a  given  threat. 
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4  Search  to  lock-on  time  MOE  computation 


The  estimation  of  the  duration  of  the  search  and  lock  operation  (dsi)  may  vary  depending 
on  several  factors.  More  precisely,  and  besides  x  and  P,  this  duration  estimate  dsi  will 
depend  on 

1.  the  FCR  characteristics,  such  as  the  beam  dimensions,  displacement  speed,  and  the 
search  pattern, 

2.  The  conditional  detection  probability  Pd\l  which  is  also  a  characteristic  of  the  FCR. 
This  measures  the  probability  that  the  FCR  detects  the  target,  given  that  the  latter 
is  within  its  gate. 

3.  The  localization  probability  distribution  p x,  which  is  defined  by  the  2D  track  infor¬ 
mation  provided  by  the  tracking  system. 

4.  The  motion  prediction  model.  Since  the  FCR  is  searching  for  non-stationary  target, 
the  2D  track  information  must  be  time-updated  using  a  dynamic  model  of  the  target 
motion. 

5.  Some  empirical  parameters  (defined  in  the  following  sections). 

All  these  factors  must  be  known  or  determined  in  order  to  estimate  the  expected  duration 
of  the  search  and  lock-on  operation.  Within  this  process  estimation,  the  possibility  of 
sweeping  all  the  search  area  without  finding  the  target  must  also  be  considered.  In  this 
case,  the  search  must  be  restarted  within  a  redefined  area  wider  than  the  previous  one. 
The  area  redefinition  must  take  into  consideration  the  possibility  that  the  target  could 
have  moved  meanwhile  and  could  be  outside  of  the  previously  defined  initial  search  area. 
Several  iterations  may  be  required  before  detecting  the  target.  In  short,  the  problem  we 
are  investigating  addresses  the  following  issues 

1.  Delimiting  the  search  area  for  the  FCR  based  on  the  information  provided  by  the 
tracking  system  associated  with  search  radars. 

2.  Determining  the  detection  and  localization  probabilities  based  on  covariance  informa¬ 
tion  provided  by  the  same  tracking  system. 

3.  Defining  a  method  for  the  estimation  of  the  search  and  lock  duration  based  on  the 
above- given  information  and  the  characteristics  of  the  underlying  FRC. 

4.  Updating  periodically  both  the  search  area  and  localization  probabilities.  This  is 
required  to  take  into  account  any  possible  movement  of  the  target  during  the  search 
operations.  This  boils  down  to  updating  the  target  state  x  and  the  covariance  matrix 
P  and  perform  tasks  1  &  2  on  the  newly  obtained  information. 

Each  of  these  issues  raises  several  sub-problems  that  will  be  discussed  in  the  next  sections. 
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4.1  FCR  search  area 


Since  the  tracking  system  associated  with  the  search  radars  provides  a  state  x  and  the 
covariance  matrix  P,  this  information  is  converted  into  probability  assuming  normal  distri¬ 
butions.  The  search  area  A  is  determined  according  to  the  target  distribution  in  bearing  (3 
along  with  a  confidence  level  a.  Since  the  2D  tracking  system  provides  no  elevation  informa¬ 
tion,  the  target  localization  probability  is  unknown  in  elevation  (e).  A  uniform  probability 
distribution  is  then  considered  with  an  elevation  range  [0°,90°].  Note  that  some  identi¬ 
fication/classification  information  may  be  used  to  reduce  this  range  based  on  the  threat 
characteristics  (e.g.,  high  diver,  sea-skimmer,  etc..). 

The  search  area  A  is  modeled  such  that  the  area  is  partitioned  into  squared  cells  having 
dimensions  equal  to  the  radar  beam  dimensions  bw  x  bw.  Each  cell  (i,j)  has  its  associ¬ 
ated  localization  probability  of  the  target  PL,(hj)-  This  localization  probability  is  defined 
according  to  the  track  information  provided  to  the  FCR  by  the  tracking  system. 

The  search  area  is  therefore  bounded  in  bearing  according  to  a  confidence  level  cq  which 
determines  the  minimum  bearing  /?_  and  the  maximum  bearing  /?+  in  which  to  search.  The 
bounds  £_  and  £+  in  elevation  are  determined  by  some  predefined  criteria  or  limitations. 
Ideally,  A  has  to  be  bounded  in  bearing  so  that  the  locations  with  a  very  low  localization 
probability  pr  are  kept  outside  the  search  area  A.  Figure  7  shows  an  example  of  an  area 
scanned  by  the  FCR. 


search  pattern 


P-  3+ 

(b)  Bearing/elevation 
plane 


Figure  7:  Area  scanned  by  the  FCR 


Let  N(j3,i 7p)  be  the  bearing  normal  distribution  (see  Section  4.2.1).  N(j3,crp)  can  be  con- 
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verted  into  a  standard  normal  distribution  N( 0, 1)  (z  scores)  using 


z  — 


P~P 


(7) 


From  the  confidence  level  <a,  one  must  find  a  lower  and  upper  bounds  that  cut  out  1-a  of 
the  area  under  the  curve.  Looking  up  in  a  z  table,  those  bounds  correspond  to  the  z  scores 
—z^ and  (see  Figure  8).  /?_  and  (3+  are  defined  according  to  —  z^i^a^  and  z^i-ay 

as  follows 


/?_=/?-  2(1^2)  <5-^  (8) 

P+=  P  +  Z(i=<*)07 3  (9) 

The  width  w  and  the  height  h  of  A  are  then  defined  by 

W  =  j3+-  13-  =  2^1^073  (10) 

h  =  £+  —  £-  (11) 


Figure  8:  Bounds  in  the  standard  normal  distribution  with  a  confidence  level  a 

4.2  Detection  and  localization  probabilities 

To  each  cell  (z,  j)  of  the  search  area  A  can  then  be  associated  a  detection  probability 
PD(hj)-  The  detection  probability  po  is  the  probability  that  the  FCR  detects  the  target 
when  aiming  at  a  particular  location  (cell).  This  probability  is  defined  as  follows: 

Pd  =Pl  ><Pd\l  (12) 

where  pl  is  the  localization  probability  and  Pd\l  !he  conditional  detection  probability  of 
the  used  sensor.  The  probability  PL,(hj)  depends  on  the  target  distribution.  The  localiza¬ 
tion  probabilities  pl  are  normalized,  so  that  their  summation  in  A  amounts  to  1.0.  The 
probability  Pd\l  depends  on  the  properties  and  performance  of  the  used  sensor  (i.e.,  the 
FCR) .  This  is  a  conditional  detection  probability  that  gives  the  probability  of  detecting  the 
target,  assuming  its  location  is  known  (pl  —  1.0). 
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4.2.1  Target  localization  probability 


As  mentioned  previously,  the  localization  probability  at  a  particular  position  is  given 
using  the  probability  density  function  of  the  target  localization.  Let  us  suppose  a  Gaussian 
target  distribution  in  bearing  (3  and  a  uniform  distribution  in  elevation  e.  The  probability 
density  function  <&(/?)  of  the  target  is  univariate  and  given  by  the  normal  distribution 

N0,vp) 


O>(P)  =  N0,ap)  = 


a/27T<7 3 


exp 


(P  ~  P) 


(13) 


In  practice,  the  support  of  the  probability  density  function  is  finite  and  must  be  bounded. 


4.2.2  Bounded-support  probability 


The  bounded  support  and  normalized  probability  density  function  $5  is  given  by 


MP)  =  { 


HP) 

Z  0HPr 

0, 


P-  <  P  <  P+ 

otherwise 


(14) 


Hence,  the  localization  probability  distribution  can  be  represented  by  the  new  defined  prob¬ 
ability  density  function 

PL  =  $6  (15) 


4.2.3  Calculating  $(/3)  from  $(x) 

In  Subsection  4.2.1,  ap  was  required  to  calculate  <&(/?).  Appendix  B  gives  all  the  required 
equations  to  calculate  the  different  conversion  of  Cartesian  information  to  polar/spherical 
information  and  vice  versa.  Also  is  given  the  equations  that  allow  for  computing  the  asso¬ 
ciated  variances. 

4.3  Average  search  time  estimation 

In  this  section,  the  expected  duration  of  the  search  and  lock  task  is  estimated.  First,  the 
case  of  single-area  search  is  considered  in  Subsection  4.3.1.  The  case  of  recursive,  multi-area 
search  is  considered  in  Subsection  4.3.2 

4.3.1  One  area  estimation 

Here  the  average  time  to  find  the  target  in  a  single  search  area  A  is  estimated.  The  general 
case  will  first  be  described,  then  the  specific  example  of  a  standard  Fixed  Swath  search 
path  with  a  rectangular  search  area  is  presented. 

4.3. 1.1  General  case 

The  average  search  time  t  in  an  area  A  is  calculated  by  summing  the  detection  probability 
for  each  cell  (i,  j)  multiplied  by  the  elapsed  time  tpcR^hj)  when  the  beam  of  the  FCR 
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goes  over  cell  (i,  j).  The  elapsed  time  depends  on  the  beam  displacement  speed  and  on  its 
pathway. 

f=EE  PD(hj)  X  tpCRihj)  (16) 

i  3 

Note  that  the  position  of  a  cell  (i,  j)  is  given  in  the  Spherical  coordinates  by 

(/?,£)  =  (/?_  +ibw  +  (17) 


4.3.1 .2  Fixed  swath  path 


Lets  consider  the  particular  example  of  the  standard  Fixed  Swath  search  path,  shown  in 
Figure  7.  The  area  to  be  searched  is  bounded  as  described  in  Section  4.2.  The  average  time 
t  of  the  standard  Fixed  Swath  search  path  is  calculated  by  adding  the  average  time  for  the 
displacements  in  bearing  and  elevation,  as  follows 


t  —  t(3  +  t£ 


(18) 


with  tp  and  t£  being  the  average  time  spent  displacing  along  in  bearing  and  elevation 
respectively.  Time  t£  is  calculated  according  to 


te 


7T - 1 

bw 

E 

3=  0 


3  r 


E 

k= 0  L 


1 


]4  Po(i,k) 
i= 0 


(19) 


where  bw/b£  is  the  beam  displacement  time  of  a  bw  angle  in  elevation.  The  time  tp  is 
calculated  by  summing  the  average  time  for  the  left  to  right  sweeps  (tir)  and  the  average 
time  for  the  right  to  left  sweeps  (tri) 


tp  —  tir  tri 


with 


where 


o<  h~bw  jv _ -I 

3  2 bw  bw 


tir 


2 


j= 0  2=0 


w  ^  . 

7~2J  + 1 
bw 


'ip  L 


w 


—  k  +  i 


is  the  cumulated  time  according  to  cell  (z,  fc),  and  with 


•  ^  h  —  bw  w  _ -i 

2  bw  bw 


tri  =  2  2J  +  1)T 


j= 0  2=0  P  *- 

where  the  cumulated  time  according  to  cell  (*,  k)  is 

b_l 

bp  L 


w  „ .  .  w 

—  {2j +  !)  +  {— -1) -i 

Uin  Oil, 


w  1  w  .. 

—k  +  {—  -1)-1 

Uin  Uii, 


(20) 

(21) 


(22) 


(23) 


(24) 
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4.3.1. 3  Gaussian  PL 


Having  a  Gaussian  localization  probability  distributions,  the  average  search  time  is  related 
to  the  target  estimated  bearing  position  j3  in  the  search  area.  The  average  time  spent 
searching  the  target  will  be  the  time  the  FCR  takes  up  to  the  estimated  position  in  the 
bearing.  With  Gaussian  distributions,  the  target  estimated  bearing  position  (/3)  is  at  the 
center  of  the  search  area.  Then,  the  sweeping  time  up  to  the  estimated  position  is  half  the 
time  to  sweep  the  area  totally.  Let  T  be  the  time  to  sweep  the  area  totally 

T={yj-K)_h  |  h  -  K  ^ 

b@  bw 


where 


{w  -  bw)  h 
bp  bw 


is  the  time  spent  displacing  along  the  bearing  axis,  while 


h  bw 
be 

is  the  time  spent  displacing  along  the  elevation  axis.  Thus,  the  average  search  time  in  one 
area  with  a  Gaussian  localization  probability  distribution  is 


£ _  T  _  (w  b'uj )  h  ^  h  bw 

2  2  bp  bw  2b£ 


(26) 


4.3.2  Recursive  estimation 


The  track  information  x  and  P  provided  to  the  FCR  by  the  tracking  system  is  used  to 
define  the  initial  search  area  A  and  the  associated  localization  probability  $5.  These  two 
quantities  allow  for  computing  the  average  time  (t)  required  to  re-acquire  the  target  by  the 
FCR,  assuming  that  the  probability  to  find  the  target  in  A  is  equal  to  1.0. 

Lets  An  be  the  area  the  FCR  has  to  sweep  at  the  77th  scan.  This  means  that  the  FCR  was 
unsuccessful  in  re-acquiring  the  target  in  its  n  —  1  previous  scans.  The  time  tn  defined  in 
Equation  16  is  then  the  average  search  time  for  scan  n  only.  Tn  is  the  time  required  to 
scan  area  An  totally,  as  defined  by  Equation  25.  At  scan  n,  if  one  is  sure  the  target  will  be 
detected,  the  required  search  time  would  be  at  least 

Tn- 1  +  Tn- 2  +  ...  +  Tf 


and,  at  most, 

The  expected  time  would  be 


Tn  +  Tn- 1  +  Tn- 2  +  ...  +  T\ 

in  +  Tn- 1  +  Tn- 2  +  ...  +  T\ 


When  the  FCR  is  cued,  one  is  able  to  calculate  Ti,  the  time  to  scan  the  initial  area  entirely 
according  to  x\  and  Pi,  the  state  estimate  and  its  corresponding  error  covariance  matrix 
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just  before  the  FCR  begins  scanning.  Ti,  x\  and  P\  then  allow  for  computing  X2  and  P25 
the  state  estimate  and  its  corresponding  error  covariance  matrix  just  before  the  FCR  begins 
the  second  scan,  using  the  state  prediction  model  (Appendix  A).  Recursively,  X2  and  P2 
will  engender  T2  and  £2,  and  so  on.  Now,  let 

1.  P]jn  be  the  probability  to  detect  the  target  at  scan  n  in  a  search  area  An,  as  given 
by  Equation  27, 

2.  Psn  be  the  probability  that  the  search  goes  through  n  trials  before  finding  the  target, 
as  defined  in  Equation  28,  and 

3.  tsn  be  the  estimated  elapsed  time  going  up  to  n  trials  to  find  the  target,  as  given  by 
Equation  29. 

Table  2  shows  the  estimated  search  times  and  the  probabilities  of  going  through  some 
number  of  trials  before  finding  the  target. 


#  Scans 

Estimated  time 

Probability 

i 

isx  ti 

pSi  =  pDi 

2 

is2  —  (^i  +  ^2) 

Ps2  =  (1  -  Pd1)Pd2 

3 

ts3  —  (Ti  +  I2  +  £3) 

pSs  =  (!  -  -Pdi)(1  -  Pd2)Pd3 

4 

ts4  —  (7 1  +  T2  +  T3  +  £4) 

psA  =  (!  -  pD i)(1  -  Pd2){  1  -  Pd3)Pda 

Table  2:  Estimated  time  and  probability  related  to  the  number  of  trials 


pDn  =  1  -  1  -  PD(i,j ) 

(i,j)eAn 
n—1 

Psn  =  Pdu  ■  JJ(1  -  Pnk ) 

k= 1 
n—1 

tsn  =tn  +  ^  Tk 
k= 1 

where  Yl(ij)eAn  stands  for  “product  on  all  cells  in  Anv .  Adding  the  probabilities  Psn  for 
each  event  n  (number  of  trials)  multiplied  by  the  related  time  spent  tsn ,  one  obtains  the 
expected  search  time  dsi  for  a  target  considering  N  trials 

N 

dsi  =  y^JsnPsn  (30) 

n=l 

Remark 

To  obtain  a  plausible  estimation  from  Equation  30,  a  number  of  iterations  N  must  be 
predetermined  or  bounded  beforehand.  To  avoid  that  probability  detection  P b  vanishes, 
the  re-acquisition  should  happen  within  few  iterations.  Otherwise,  the  search  areas  will 


(27) 

(28) 
(29) 
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grow  unboundedly  and  become  too  large  (i.e.,  the  localization  probability,  too  low)  to 
detect  any  object. 

Although  tsn  (and  the  search  area  A)  grows  exponentially  and  can  become  intractable 
very  fast  (exponentially)  for  some  particular  initial  parameters,  Psn  will  vanish  rapidly,  so 
that  the  solution  should  converge  as  the  number  of  trial  increases.  The  periods  of  time 
corresponding  to  the  first  trials  should  be  the  most  contributive  to  the  solution. 

4.3.3  Stop  conditions 

As  mentioned  previously,  the  process  of  search  and  lock-on  may  involve  many  sweeps  in 
order  to  re- acquire  the  target  and  lock  on  it.  The  process  could  even  go  non-stop  and 
limitless  if  the  FCR  never  re- acquires  the  target.  For  the  purpose  of  time  estimation,  the 
process  must  be  stopped  at  some  instant  that  can  be  defined  as  a  specific  number  of  sweeps 
(or  iterations).  The  problem  then  comes  to  find  how  the  number  of  sweeps  will  be  set.  The 
solution  will  depend  on  the  application,  such  as  those  covered  with  the  problems  mentioned 
in  Subsection  3.3.1. 

However,  besides  the  application  related  issues,  other  factors  may  set  some  limits  to  the 
number  of  sweeps.  One  factor  is  the  workload  involved.  In  the  search  process,  the  search 
area  grows  exponentially  at  each  trial  and  may  become  very  large  after  some  number  of 
sweeps.  It  is  possible  that  the  search  area  becomes  so  wide  that  the  estimated  search  to 
lock-on  time  begins  to  grow  with  the  number  of  iterations  considered  instead  of  stabilizing 
to  a  specific  value.  This  may  lead  to  an  infinite  estimated  search  to  lock-on  time. 

In  such  a  situation,  t$n  is  so  large  that,  even  though  Psn  becomes  smaller  as  the  number  of 
iterations  (n)  increases,  the  product  tsnPsn  (in  Equation  30)  increases  instead  of  decreasing 
towards  zero.  When  such  a  situation  is  detected,  the  iterative  estimation  should  be  stopped 
to  avoid  unnecessary  calculus. 

As  an  example  of  a  stop  condition,  one  can  avoid  excessive  calculations  by  setting  the  total 
detection  probability  required  PDtot  •  The  total  detection  probability  PDtot  is  related  to  the 
number  of  trials  N  considered  (Equation  31).  Having  PDtot  =  1  would  require  an  infinite 
number  of  trials.  Thus,  by  setting  PDtot  <  1?  one  gets  a  finite  number  of  trials  N 


N 


PDtot  ~  1  Pmtot  ~  1  JJ  Prr 


N 


»-n 


n  0--pd{}J)) 

(ij)eAn 


(31) 


where  Pmtot  is  the  total  miss  probability  and  Pmn  is  the  miss  probability  at  trial  n.  Note  that 
PD^hj)  is  not  the  same  in  each  trial  n.  In  practice,  Equation  31  is  calculated  iteratively  by 
considering  each  trial  n,  and  thereby  the  iteration  should  stop  when  Pr>tot  reaches  a  specific 
value. 
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5  Conclusion 


The  Decision  Support  Systems  (DSS)  Section  at  Defence  Research  &  Development  Canada 
(DRDC)  -  Valcartier  is  conducting  research  activities  aiming  at  defining,  developing,  and 
demonstrating  advanced  data  fusion  adaptation  concepts.  The  objective  of  the  reported 
work,  as  part  of  these  research  activities,  is  to  propose  a  new  Measure  Of  Effectiveness 
(MOE)  based  on  the  search  to  lock-on  time  of  the  Fire  Control  Radar  (FCR).  This  proposed 
Measure  Of  Effectiveness  aims  at  providing  the  required  adaptation  context  of  the  target 
tracking  process  in  naval  Above  Water  Warfare.  This  MOE  establishes  a  quantitative 
relationship  between  the  maritime  tactical  picture  quality  and  the  search  to  lock-on  time 
of  the  FCR. 

A  method  to  estimate  the  expected  duration  of  the  search  and  lock-on  operation  has  been 
developed.  This  estimation  takes  into  consideration  the  FCR  search  pattern,  the  beam 
width,  and  the  displacement  speed.  It  also  considers  the  localization  and  detection  prob¬ 
abilities.  The  method  iterates  over  a  series  of  sweeps  the  FCR  can  go  through  before 
re- acquiring  and  locking  on  the  target.  These  sweeps  must  be  limited  following  specific 
conditions  such  as  the  total  detection  probability  cumulated  after  each  sweep. 

Within  the  framework  of  the  adaptive  tracking  system  of  the  FCR,  the  proposed  MOE 
being  used  tackles  and  illustrates  different  adaptation  problems.  This  represents  a  first 
step  towards  building  a  tracking  system  that  meets  all  the  adaptive  requirements. 


DRDC  Valcartier  TR  2003-284 


21 


References 


[1]  Roy,  J.  and  Bosse,  E  (1995),  Definition  of  a  Performance  Evaluation  Methology  for 
Sensor  Data  Fusion  Systems.  DREV  TR  9423.  UNCLASSIFIED. 

[2]  Benaskeur,  Abder  Rezak  and  Bosse,  Eloi  (2003),  Adaptive  Data  Fusion  Concepts  for 
Tracking  in  a  Tactical  Situation,  In  RTO  -  SET  059  Symposium  on  Target  Tracking  & 
Data  Fusion  for  Military  Observation  Systems. 


22 


DRDC  Valcartier  TR  2003-284 


Annex  A:  Prediction  model 


The  tracked  target  is  assumed  to  be  moving  in  a  2D  space,  where  the  acceleration  acts  as 
an  input5.  Cartesian  coordinates  are  used  in  order  to  be  in  a  linear  system.  The  state  to  be 
estimated  is  composed  of  the  target’s  coordinates,  viz .,  the  position  and  the  linear  velocity. 
With  the  following  state  variable  notation 


x  — 


and 


v  = 


the  equations  of  such  a  target  can  be  expressed  as 


(A.l) 


xk+l  —  F'kxk  “1“ 


(A- 2) 


where  v  is  vector  of  random  variables  that  reflects  the  unforeseeable  variation  of  the  accel¬ 
eration  in  both  directions  and  the  state  transition  matrix  is  given  by 


Fk 


1  0  At  O' 
0  1  0  At 

0  0  1  0 
0  0  0  1 


(A.3) 


where  At  is  the  time  increment.  The  matrix  T  in  (A. 2)  depends  on  the  model  used  to 
represent  the  discrete-time  nature  of  the  process  noise  v.  Two  examples  are  given  below. 
The  first  one 


"At2/2  0 

0  At2/2 

At  0 

0  At 


(A.4) 


represents  to  the  pulse  model,  while  the  following  one 


ft/At2 


r  = 


2VS 

0 


0 

0 


v^At2 

2 

0 

\/A t 
0 


0  0 

v^At2  ^At2 

2y/S  2 

0  0 

0  v/A t 


(A.5) 


is  used  by  the  Brownian  motion  modeling.  Note  that  the  latter  requires  four  random 
variables  instead  of  two.  If  the  process  noises  are  assumed  to  have  the  same  standard 

5 Since  the  acceleration  may  change  unforeseeably,  it  is  often  modeled  as  a  random  variable,  and  so  will 
be  as  such. 
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deviation  dv ,  the  pulse  model  results  in  the  following  process  noise  covariance  matrix 


Q  =  S2vTTt 

'At4 


=  St 


4 

0 

At3 


2 

0 


0 

At4 

0 

At3 

~Y 


At3 

~Y 

0 

At2 

0 


0 

At3 

— 

0 

At2 


while  the  covariance  matrix  for  the  brownian  motion  is  given  by 


Q 


<52rrT 

'At3 


3 

0 

At2 

0 


0 

At3 

0 

At2 

~Y 


At 2 

~Y 

o 

At 

0 


0 

At2 

~Y 

o 

At 


(A. 6) 


(A. 7) 


(A.8) 


(A.9) 


We  assume  that  a  Kalman  filter  is  used  for  tracking.  Since  there  is  no  sensor  report  during 
the  search  and  lock-on  phase,  only  the  time  update  part  of  the  filter  is  considered.  The 
state  dynamic  model  (Equation  A. 2)  serves  for  the  prediction  of  the  state  variable  x  and 
the  corresponding  covariance  matrix  at  trial  n  +  1  from  their  values  at  trial  n  and  Tn. 

Note  that  Tn  is  the  time  to  sweep  the  area  An  totally  at  scan  n  (Equation  25). 

1.  Prediction  of  the  state  —  The  most  recent  available  a  posteriori  estimate  of  the 
state,  namely  xTn\Tn ,  is  injected  in  the  system  model  to  yield  the  new  a  priori  estimate 
xk+1\k  of  the  state  at  the  time  instant  Tn+ 

*Tn+i|T„  =  FTnXTn\ Tn  (A.  10) 

2.  Prediction  of  the  covariance  matrix  —  A  prediction  of  the  estimation  error 
covariance  matrix  is  performed,  at  the  time  instance  Tn+ Its  computation  is  based 
upon  the  most  recent  available  a  posteriori  value  PTn\Tn. 

PTn+i\Tn  =  FT„PTn\TnFTn  +  QT 

Note  that  the  system  noise  v  is  not  used  in  the  state  prediction  equation  (A.  10).  This  is 
due  to  the  fact  that  the  best  prediction  one  can  make  for  a  white  noise  is  its  mean,  which 
happens  to  be  equal  to  zero  here. 
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Annex  B:  Coordinate  systems 


The  spherical  coordinate  system  family  allows  for  tracking  in  the  same  coordinate  system 
as  the  sensor  reports.  It  is  particularly  more  suitable  when  the  sensor  reports  only  partial 
information  (e.g.,  angle  only).  For  the  widely  used  constant  velocity  and  constant  accel¬ 
eration  models  of  the  target  flight  trajectories,  tracking  in  the  spherical  system  makes  the 
motion  equations  lead  to  highly  nonlinear  systems.  In  the  following,  when  required,  the 
spherical  coordinate  system  (i?,  0,  </>)  is  replaced  by  a  slightly  modified  version,  hereafter 
simply  referred  to  as  the  spherical  coordinate  system  (i?,  /?,  e).  R  is  a  distance  coordinate 
that  corresponds  to  the  slant  range  relative  to  the  interceptor,  (3  (bearing)  is  defined,  in 
the  horizontal  plan,  as  the  angle  relative  to  the  North,  while  e  (elevation)  is  defined  as  the 
elevation  relative  to  the  Horizon  (see  Figure  B.l).  This  modification  is  only  introduced  to 
conserve  the  compatibility  with  the  previous  works. 


Up  (z) 


Figure  B.1:  Spherical  coordinate  system  (adapted  convention) 


The  Cartesian  coordinate  system,  on  the  other  hand,  represents  the  natural  choice  for  the 
widely  used  constant  velocity  and  constant  acceleration  dynamical  models  of  the  target 
flight  trajectories.  The  main  advantage  of  this  system  lies  in  that  it  allows  for  the  use 
of  the  well-known  standard  Kalman  filtering  technique,  which  is  applicable  only  when  the 
state  transition  and  observation  equations  are  expressed  by  linear  functions. 
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B.1  Cartesian  3D  to  spherical  conversion 


Based  on  the  available  Cartesian  information,  the  spherical  coordinates  and  their  derivatives 
are  expressed  as 

R  =  {x2  +  y2  +  z2)1/2  (B.l) 

R  =  (xx  +  yy  +  zz)(x2  +  y2  +  z2)~1^2  (B.2) 

(3  =  tan-1  xy -1  (B.3) 

P  =  (yx  -  xy)(x2  +  y2)-1  (B.4) 

e  =  tan-1  z(x2  +  y2)-1/2  (B.5) 

e  =  ( x 2  +  y2)-1/2  z  —  z{xx  +  yy  +  zz)(x2  +  y2  +  z2)-1  (B.6) 


Each  of  the  above-given  variables  can  be  rewritten  as 


dA 

dt,  d£ 

dt, 

9£l 

_ _  y 

_dx 

dy 

dz  dx 

dy 

dz\ 

X  C 


The  corresponding  error  variances/covariances  are  given  by 


=  E\K1  —  Jt  E  XX 1  L7C 


"2  2  2 
®x  ®xy  ®xz 

a2  a2 


az  ®zx  azy  azz 

<A  at-,  ah 


ay  ayz 


which  can  be  rewritten,  using  the  partial  derivatives,  as 

_2  d£dC_2  ,  d^d(  2,d£dC_2  ,  dtd(_2 


2  _  di 


2  , 

9^*  + 
d£d(_ 
dx  dy 

<9£  <9C 

dx  dy 
d^dC 
dy  dx 

dz  dx 
d^dC 
dx  dy 


__2  | 
axy  + 


2  , 
®yx  ' 


azx  + 


9a;  9y 

tKdt 

dx  dy 

(K(K 

dy  dx 
<9C<9£ 
dz  dx 

tKdt 

dx  dy 


2  i 

axy  + 


5C  2  ,  ; 

dz  °z  dx  dx  G' 
d. £  d(  9£  dl ; 

dx  dz  dx  dz 

9£  9£  9£  9£ 

9a;  9i  9a;  dz 

'9C9C  9C9e 

dy  dy  dy  dy 

'dt[(X 

dz  dy  dz  dy 

9£  9£  9£  9£ 

dx  dz  dx  dz 


axz  + 


2  I 

“T 


+  2  , 

9y  9y  ^ 

2  ,  \dtdC 

°xz  [dx  dx 

2  , 

CTa:i  dydz 

2  , 

ayy  dydz 

2  , 

[dzdz 

2  , 

xz  [dydz 


axz  + 


0£0C  2 

dz  dAz 

JXd£ 

dx  dx 

dc^ 

dy  dz 
dCdt~ 
dy  dz 
9C9^: 
dz  dz 

.dcoe 

dy  dz 
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For  each  spherical  dimension,  the  different  partial  derivatives  are  given  below,  where  we 
define 


R2D  =  ( X 2  +  y5 


'dR 

dR 

dR 

dR 

dR 

dR 

dx  ’ 

v 

dz 

’  d±' 

w 

dz 

r 

1 T 

xR~ 

\  yR 

-l 

•> 

zR~\ 

0,  0, 

,  0 

'dp 

dp 

d/3 

dp 

dp 

dp V 

dx  ’ 

v 

dz 

’  d±' 

w 

dz 

yRJ>  ~xRon’ 


0,  0,  0,  0 


(B.10) 

(B.ll) 


\de 

de 

de 

de 

de 

del 

dx 5 

dy ' 

dz' 

dx 

w 

dz 

-xzR^R  2,  -yzR2*R  2,  R2DR  2,  0,  0,  0 


(B.12) 

(B.13) 


dR 

di? 

dR 

dR 

dR 

dR 

dx ’ 

ay’ 

d? 

dx 

dy 

dz 

xR 2 

—  X\ 

[xx  + 

yy  + 

zz) 

y-R2 

-  y(xx  + 

yy  + 

zz) 

ii?2 

-*( 

'xi;  + 

yy  + 

zz) 

xR 2 

XL 

yR2 

zR2 

(B- 14) 


(B.15) 


dp  dp  dp  dp  dp  dplT 
dx  dy  dz ’  dx  dy  dz 
-yR22D  ~  2 x(xy  -  xy)~ 
xR\d  ~  2 y(xy  -  xy) 


(B.16) 


(B.  17) 
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(B.18) 


Jc  = 


"  de 

de 

de 

de 

de 

de 

T 

_dx 

3/ 

dz 

dx 

3/ 

dz_ 

xz(xx  +  yy)Rj  +  (xz  —  xz)  —  2xzR ; 


-2D  ■  V—  - 2DR  2  +  2xz(xx  +  yy)R 

yz(xx  +  yy)R~p  +  ( yz  -  yz)  -  2yzR2D 

»2  D  —  2  I  o  ~2 


2D 

-(xx  +  yy) 


2zzRz2DR~ 

—xz 

-yz 

RId 


~2  +  2  xz(xx  +  yy)R~ 21 
R~2  +  2  yz(xx  +  yy)R~2 
+  2  z2(xx  +  yy)R~2 


R 


—2  p— 1 
' 2£> 


(B.19) 


B.2  Spherical  to  cartesian  3D  conversion 

In  the  way,  based  on  the  assumed  available  spherical  information,  the  Cartesian  coordinates 
and  their  derivatives  are  expressed  as  follows 


x  —  R  sin  f3  cos  e 
y  —  R  cos  / 3  cos  e 
z  —  R  sin  e 

x  —  R  sin  (3  cos  e  +  R(3  cos  (3  cos  e  —  Re  sin  / 3  sin  e 
y  —  R  cos  f3  cos  e  —  i?/?  sin  (3  cos  e  —  Re  cos  /?  sin  e 
z  —  R  sin  e  +  Re  cos  e 


and  rewritten  in  a  compact  form  as 


\dd 

dd 

dd 

dd 

dd 

dd] 

OR 

dp 

de 

dR 

dp 

de 

~R~ 

p 


which  leads  the  following  error  variances/covariances  expression 
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/ 3R 
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a 


a 


a 
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3c, 


(B.20) 

(B.21) 

(B.22) 

(B.23) 

(B.24) 

(B.25) 


(B.26) 
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This  can  be  rewritten  using  the  partial  derivatives  as 
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where  the  different  partial  derivatives  are  given  below. 


Jx  = 


dx  dx  dx  dx  dx  dx 

OR’  dp '  Se’  Cr’  d/C  Si 


sin  P  cos  e 
R cos  P cos  e 
— R  sin/?  sine 
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(B.27) 
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dy 
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— R  cos  (3  sine 
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de 

0 
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0 
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(B.28) 


Jz  = 
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(B.30) 


Jx  = 
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(B.36) 

(B.37) 


(B.38) 


B.3  Cartesian  2D  to  polar 

The  polar  coordinate  are  a  2D  special  case  of  the  spherical  coordinates.  The  coodinate 
transformations  from  Cartesian  (x,y)  to  polar  (R,  ft)  are  derived  from  the  spherical  as 
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follows 


R  —  \J  x1  +  y2 

[3  —  tan-1  — 

V 

xx  +  yy 

\/?+? 
yx  —  xy 


R  = 


ft  9  i  9 

xz  +  yz 

The  above  given  variables  can  be  rewritten  as 


£  = 


and,  as  in  the  previous  case, 
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which  can  be  rewritten  using  the  partial  derivatives  as 
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where  the  partial  derivatives,  of  each  dimension,  are  given  below. 
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(B.44) 

(B.45) 
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(B.47) 
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Ju  = 


dR  OR  dR  dk\T 
dx'  dy ’  dx  dy 
xR 2  —  x2x  —  xyy 
ijR2  -  xyx  -  y2y 

xR2  R 
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(B.48) 


(B.49) 


dp  dp  dp  dplT 
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yR 2 

—xR2 


(B.50) 


(B.51) 


B.4  Polar  to  cartesian  2D 

In  the  2D  case,  the  Cartesian  coordinates  and  their  derivatives  are  obtained  as  follows: 


x  —  R  sin  (3 
y  —  R  cos  P 
x  —  k  sin  P  +  pR  cos  P 
y  =  k  cos  P  —  PR  sin  P 

Each  of  the  above  given  variables  can  be  rewritten  as 


(B.52) 

(B.53) 

(B.54) 

(B.55) 


OS,  dt,  d  g  dp 

dR  dp  dR  dp. 


=  jjx 


(B.56) 


with  as  a  variance/covariance  matrix 


‘2  —  ST  l 


whose  development  is  given  by 
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dR  dR°R  +  dp  dpa/3  +  dR  dR °R  +  dp  dp° 
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where  the  Jacobians  are  given  by 


Jx  —  [sin/3,  i?cos/3,  0,  0]T  (B.57) 

Jy  —  [cos/3,  — i?sin/3,  0,  0]T  (B.58) 

Jx—[p  cos/3,  Rcos/3  —  /3i?sin/3,  sin/3,  Rcos/3]T  (B.59) 

=  [ — /?  sin  /3,  — i?sin/3  — /3i?cos/3,  cos/3,  —  Rsm/3]T  (B.60) 


B.5  Polar  contact  to  cartesian 


In  the  case  of  polar  contacts,  the  only  available  information  is  the  range  R  and  the  beating 
/3,  i.e.,  no  information  about  the  rates  is  available.  Hence  the  possible  transformations  are 


x  —  R  sin  /3 
y  —  R  cos  /3 


This  reduces  the  above-given  notation  to 

with  as  a  variance/covariance  matrix 
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dRdRaR  +  df5dpap 


which  gives  for  x  and  y 


a2  =  sin2  /3a2R  +  R2  cos2  /3a p 
a2  =  cos2  fla  ji  +  R2  sin2  [3 a2 
a2y  =  sin  / 3  cos  P(aR  —  -R2<t|) 


(B.61) 

(B.62) 


(B.63) 


(B.64) 

(B.65) 

(B.66) 
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